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We present a theoretical study of the energetics of thin nematic shells with two charge one-half 
defects and one charge-one defect. We determine the optimal arrangement: the defects are located 
on a great circle at the vertices of an isosceles triangle with angles of 66° at the charge one-half 
defects and a distinct angle of 48°, consistent with experimental findings. We also analyse thermal 
fluctuations around this ground state and estimate the energy as a function of thickness. We find 
that the energy of the three-defect shell is close to the energy of other known configurations having 
two charge-one and four charge one-half defects. This finding, together with the large energy barriers 
separating one configuration from the others, explains their observation in experiments as well as 
theirlong-time stability. 


I. INTRODUCTION 

One of todays major drives in condensed matter 
physics is the assembly of mesoscale particles into com¬ 
plex structures [T]. By creating anisotropy in the in¬ 
terparticle interactions, one can increase the complexity 
and functionality of these structures. A proposed way to 
achieve anisotropic interactions is by coating a spherical 
particle or droplet with an orientationally ordered phase 
[2]. The topology of the sphere enforces defects in the 
coating. Since these defects are very distinct regions on 
the sphere, they are suitable for the attachment of link¬ 
ers acting as bonds between the particles. For the case of 
a vector order parameter, topology requires two defects, 
creating a particle with two binding sites. Exploiting 
this, de Vries et al. succesfully assembled chains of such 
divalent nanoparticles [3]. Nematic rather than vector or¬ 
der allows for defects of charge one-half, referring to the 
180 degrees rotation experienced by the local average ori¬ 
entation of the nematic molecules, n, when encircling the 
defect. In fact, it is energetically favourable for defects of 
charge one to split into two charge one-half defects (Fig. 
Ef)- Nematic order on the sphere has four topological de¬ 
fects of charge one half in its ground state, such that the 
sum of all charges is equal to 2, the Euler characteristic 
of the sphere, as demanded by the Poincare-Hopf the¬ 
orem. Their mutual repulsion drives them as far away 
from each other as possible: at the vertices of a regu¬ 
lar tetrahedron [4]. Thus, chemical functionalisation of 
the defects in the ground-state of two-dimensional ne¬ 
matic liquid crystal on the sphere would thus result in 
the diamond structure [2]. In the decade that followed 
the conception of this idea, a vast amount of theoretical 
and numerical work was performed; it included studying 
the effect of elastic anisotropies, external fields, spheric¬ 
ity or shape, and thickness and thickness inhomogeneity 
of the nematic film [SHIS]- These works were mainly fo¬ 
cussing on spherical nematic shells with four defects with 
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a charge of one-half. Experimental investigations on ne¬ 
matic shells generated by trapping a water droplet inside 
a nematic droplet, however, have revealed the existence of 
a much wider variety of defect structures besides the reg¬ 
ular tetrahedral defect arrangement [6l [T2j [181 EH |30H32] , 
even with a valence number different from fonrp [T2l[^ . 
There exist divalent configurations in which instead of 
four charge half defect lines spanning the shell, there are 
two pairs of point defects, called boojurns, residing on 
the boundary surfaces. They arise because the thick¬ 
ness of the nematic coating is nonzero: the elastic energy 
of a singular line with a winding number of one at the 
boundary is reduced by escaping in the third dimension, 
as illustrated in Eig. [^. This route thus forms an al¬ 
ternative to the splitting into s = 1/2 lines spanning the 
shell. 

Surprisingly, also structures containing both boo j urns 
as well as charge one-half disclination lines have been 
observed in experiments ina. These defects structures 
have threefold valence yet they are still consistent with 
Poincare-Hopf’s theorem, because the total topological 
charge of the defects at the boundary is 1 + 1/24-1/2 = 2, 
the Euler characteristic of the sphere. Again, also this 
energetically stable defect configuration arises because of 
the finite thickness of the nematic coating. If the shell is 
strongly inhomogeneous in thickness, experiments show 
that the defects cluster in the thin part of the shell mm- 
In simulations, the trivalent state has also been observed 
for inhomogeneous shells [25]. The defects were found 
to be positioned at the vertices of an isosceles triangle, 
with the boo j urns located in the thickest region of the 
shell and the s = 1/2 disclination lines in the thinnest 
hemisphere. 

In this report, we will study theoretically the defect 
separations, energetics and fidelity of the bonds in homo¬ 
geneous spherical nematic shells with three-fold valence. 
We will make a comparison with divalent and tetravalent 
shells and find the optimal valency as a function of shell 
thickness as well as the energy barrier between shells of 
different valency. We will compare these results to exper¬ 
iments in which we decrease the thickness inhomogeneity. 
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FIG. 1: (a) In a two-dimensional nematic a s = 1 
topological defect (black dot in left panel) can lower its 
elastic energy by splitting into two s = 1/2 defects 
(purple dots in right panel), (b) A singular line (left 
panel) spanning the shell with a winding number of one 
at the boundaries is topologically and energetically 
unstable. The singular core is indicated by a black dot 
in the top view shown in the top panel and by the 
vertical bold line in a cut shown in the bottom panel. 
One way of reducing the elastic energy escaping in the 
third (vertical) dimension (right panel), thereby leaving 
a point defect (green dot), called boojum, on each 
boundary. 


where Si is the topological charge of defect i, /3ij is the 
angular distance between defects i and j, and Z is the 
number of defects or valence number. The self-energy 
Ei {R) is given by 


Ei{R) = TrKsi log (E^+E„ (3) 

where R is the radius of the sphere and a is a small 
scale cut-off preventing a divergence of the energy. Ec 
represents a core energy, which depends on the details 
of the microscopic interactions. The self-energy is re¬ 
sponsible for the splitting of +1 defect in an ideal two- 
dimensional nematic, because of its proportionality with 
s‘f. The other term in eq. (H describes the repulsion be¬ 
tween like-charged defects. We wish to find the optimal 
location for the defects in a thin homogeneous shell given 
that 5i = 1, 52 = I and 53 = ^. This requires minimis¬ 
ing the interaction term of the free energy. We minimise 
the interaction energy with respect to three independent 
variables, namely / 3 i 2 , Pis and the angle, C, subtended 
by the two curved triangular sides (circular arcs) meeting 
at the charge one defect. If we apply the law of cosines 
on the sphere: 


cos P 23 = cos P 12 cos pis + sin P 12 sin Pis cos C, (4) 

we can eliminate P2S in favour of C in the free energy, 
and demand ~ ~ From the latter 

equation, = 0 , we obtain C = tt, implying that the 
defects lie on a great circle (see Figs. [^and[^. There is 
always a circle that can be drawn through three points 
on a sphere; the maximal radius of this circle reflects the 
repulsive nature of the defects. With some straightfor¬ 
ward algebra the other two equations, and = 0 , 
then lead to 


II. TRIVALENT GROUND STATE 


P12 = Pi^ = TT — arccos - « O.TStt « 131.8° (5) 

o 


The free energy of a thin curved nematic film is 


P 2 S = 2 arccos - ^ 0 . 547 r 96.4° 

o 


(6) 


F = 


^jdAki + h {DiTij - DjUi) {D^n^ - Dh 

( 1 ) 

with ki and ks the two-dimensional splay and bend elas¬ 
tic constants and Di is the covariant derivative. Eq. 0 
can be recast in terms of defect separation rather than 
the director field n. For a spherical nematic, the elastic 
energy in the one-constant approximation k = ki = ks 
reads 


z z 

F = -TTk'^SiSj log(l - cos^ij) + (R) (2) 

i<j i 


We thus find that the defects are located at the vertices 
-|Of an isosceles triangle rather than equilateral triangle, 
) shown in Figs. [^and|^ This less symmetric configura¬ 
tion arises because of the asymmetry in the magnitude 
of the charges of the defect: the two +1/2 defects repel 
each other less strongly than a charge one and charge 
one-half such that P 12 and Pis are larger than P 2 S‘ This 
is in marked contrast with the regular tetrahedral config¬ 
uration in which all the defects are equidistant, because 
all four charges are indistinguishable. The fact that 52 
and 53 are of equal magnitude is still reflected in the 
equal length of two of the sides {P 12 = P23) of triangle. 
Perhaps surprisingly, the distance between two charge 
one-half defects is smaller in the trivalent state than in 
the more ‘crowded’ tetravalent state. The surface angles 
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FIG. 2: Four views on the bend texture of the director 
field on the sphere containing a +1 defect and two +1/2 
defects arranged in an isosceles triangle with 
/3 i 2 = As « 132°, As ^ 96.4°, « 48° and 

0^2 = 0^3 ~ 66°. The defects lie on a great circle. 


FIG. 3: Four views on the splay texture of the director 
field on the sphere containing a +1 defect and two +1/2 
defects arranged in an isosceles triangle with 
A2 = As ^ 132°, As ~ 96.4°, « 48° and 

0^2 = 0^3 ~ 66°. The defects lie on a great circle. 


of the flat triangle can be found be simple trigonometry: 
by realising that the triangle formed by two defects and 
the centre of the sphere is also an isosceles triangle (of 
which two sides have a length equal to the radius) we 
obtain 


Qfi = TT — A2 = arccos - 
o 


0^2 = 0^3 = 


A2 


TT 

2 


arccos 


48.2° 
2 


: 65.9° 


( 7 ) 

(8) 


Given the defect locations the energy-minimising di¬ 
rector field can be found by means of a stereographic 
projection of the planar solution. The bend texture is 
displayed in Fig. Rotating this director field over an 
angle a yields the same free energy in the one-constant 
approximation. The splay texture (Fig. corresponds 
to a = 7r/2. We note that the escape in the third di¬ 
mension, in which the singular region is distributed over 
a larger distance of the order of the thickness, occurs in 
shells of finite thickness and is somewhat different than 
the problem of three point defects in a two-dimensional 
nematic solved above. However, we expect that the de¬ 
fect separations will be marginally affected as long as the 
thickness is small compared to the radius. 

To test the theoretical expectations, we generate ne¬ 
matic shells using microfluidics [6]. The shells produced 


with this method are double emulsions where an inner 
aqueous droplet is contained inside an outer liquid crys¬ 
tal droplet which is, in turn, dispersed in an aqueous 
solution. We use salt to establish the osmotic pressures 
of both the inner droplet and the continuous phase and 
thus the osmotic pressure difference between them. The 
stability of the emulsion is guaranteed by the presence of 
polyvinyl alcohol (PVA), which also enforces planar an¬ 
choring of the liquid crystal, 4-Gyano-4’-pentylbiphenyl 
(5GB), at both the inner and outer interfaces. The shells 
are heterogeneous in thickness due to buoyancy effects; 
they are thinner at the top and thickest at the bottom 
[6]. Hence, the shell thickness gradually increases from 
the top to the bottom of the shell. Typical values of 
the outer radius, R, and thickness of the shell, h, are 
in the ranges [20,60]/im and [1,10]/im, respectively. In 
this shell geometry, the defects are confined to the top, 
as shown in Fig. |^a), since this reduces the Frank free 
energy m- 

As a result of the imposed osmotic pressure differene 
between the inner droplet and the continuous phase, the 
inner droplet swells while the shell becomes thinner and 
more homogeneous. This happens quasi-statically [12]. 
We then monitor the evolution of the defects through¬ 
out the process. We find that the defects progressively 
spread, as shown in Figs. [^a-c). Interestingly, the shape 
of the triangle defined by the positions of the three de¬ 
fects changes shape as this happens; the surface angles 
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0 ^ 1 , ol 2 and 0^3 progressively evolve as the shell becomes 
thinner and more homogeneous, as shown in Fig. Bd). 
For shells with h/r < 0.03, the surface angles reach val¬ 
ues that are consistent with those predicted theoretically: 
ai ^ 46° and 0^2 = 0^3 ~ 68 °. This configuration corre¬ 
sponds to an isosceles triangle with the lower angle placed 
at the +1 defect. However, the defects are not yet in a 
great circle. In fact, for h/R = 0.03, the defects all lie in 
the upper hemisphere of the shell and are contained in 
a plane that is perpendicular to the gravitational direc¬ 
tion, g] the angle Oz between the normal of this plane, 
TV, and g is zero. As h/R becomes smaller than 0.03, Oz 
increases [see Fig. [^e)] and the defects progressively sep¬ 
arate from each other while maintaing the values of the 
surface angles and hence the shape of the isosceles trian¬ 
gle. As this happens, the shell progressively approaches 
the expected configuration for an infinitely thin nematic 
shell, where the defects lie on a great circle. Indeed, for 
h/R ~ 0.01, the +1 defect is located at the top of the 
shell, see FigJ^a), while the + 1/2 defects are in the lower 
hemisphere of the shell, as indicated with the arrows in 
FigJ^b) and (c). In this configuration, the distributions 
for the anbles P 12 and P 23 are both Gaussians chacterized 
by a mean of P 12 = 118° and ^23 = 77° and correspod- 
ing widths of Af 3 i 2 = 19° and A /323 = 15°, as shown in 
Figs. I^d) and (e). These values are slightly lower than 
the theoretical ones, indicating that the defects have not 
completely reached yet a great circle, possibly due to a 
remaining thickness inhomogeneity in the shell. Finally, 
note that throughout this process the pair of boojurns is 
located in the thinnest hemisphere, whereas in the sim¬ 
ulations of ref. [25] the pair of boo j urns was found to be 
located in the thickest hemisphere. These simulations, 
however, are carried out for shells with a much larger 
thickness than the shells in this experiment. 


III. VALENCE TRANSITIONS 


We will now proceed with an estimate of the energy 
of the trivalent shell when this escape is taken into ac¬ 
count. In doing so, we follow the arguments in ref. 
[ 5 ]. We first consider the energy when three singular 
lines are spanning the shell at angular distances reported 
above. We estimate this energy as the product of the 
two-dimensional result and the thickness, h, thus effec¬ 
tively taking k = Kh: 


(a) 

(b) 
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FIG. 4: (a-c) Cross-polarized images of a shell at 
different stages of a swelling process, which makes the 
shell more homogeneous in thickness. The shell has 
three topological defects: one defect of charge + 1 , 
characterized by four black brushes, and two defects of 
charge + 1 / 2 , characterized each one by two black 
brushes. The three defects depict a triangle whose 
angles, o^i, 0^2 and 0 ^ 3 , change as the shell swells, see 
evolution from left to right. The exact values of o^i, 0^2 
and 0^3 as a function of the shell average thickness h/R 
are shown in (d). The triangle depicted by the defects is 
initially oriented perpendicularly to the gravitational 
direction, but it tilts off as the shell swells, as shown in 
(d), where Oz stands for the angle between gravity g 
and the flat triangle normal TV. 


Ez =3 = irKh 


R 


1 


log ( -^ ) + - log ( — ) + 3.65 


R 


2E, 


nKh 

( 10 ) 

By comparing this to the energy of a shell with four 
disclination lines 


Ez =4 = irKh 



- 0.43 + 


4EZ 

wKh 


( 11 ) 


we can find the critical value for h above which the triva¬ 
lent defect configuration is energetically preferable over 
the tetravalent one: 


E'z=3 = T^Kh 


(^1 + 2 X 



0.54-h 


3Ee 

nKh 


( 9 ) 

A heuristic yet adequate method to include the escape 
is by replacing the microscopic cut-off by the thickness of 
the shell, since the singular core is spread out over spatial 
dimensions of the order of h. To account for the pair of 
boojums an energy A.2KTvh is added [5l 1^ l33] . We then 
obtain 


hlJR = e4-08-2So/^if/» / “ (12) 

V R 

Similarly, one can find the critical value for h below 
which the trivalent defect configuration is energeticaiiy 
preferabie over the divaient one by setting Ez =3 equai to 
the approximation of the energy of a sheii with two dia- 
metricaiiy opposite pairs of surface defects, F^z= 2 - Again, 
we wiii first find the energy of a sheii with two singuiar 
iines 
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FIG. 5: Defect configuration of a very thin and almost 
homogeneous trivalent nematic shell, (a-c) Cross 
polarised images of the shell at different focal planes: 
the +1 defect is in an upper plane shown in (a), while 
the + 1/2 defects are at lower planes, see arrows in (b) 
and (c). (d-e) Histograms of the angular distances 
between defects. The three defects depict an isosceles 
triangle where the two equal sides correspond to the 
distance between the +1 defect and each of the + 1/2 
defects, pi 2 = (118 ± 19)°, and the unequal side 
corresponds to the distance between the two + 1/2 
defects, P 23 = (77 + 15)°. 


^z=2 = '^Kh 



- 0.69 + 


2E, 

irKh 


(13) 


after which we apply the same trick as we used to find 
Ez =3 to obtain 


Ez=2 = TrKh 



+ 7.69 


We find a very similar value 


(14) 


ll;3/i? = e4•04-2BcAX/^ /« (15) 

V R 

The energy as a function of thickness is plotted in Fig. 
[ 6 ] for all three different valencies. 

Since ^23 < ^34 there is no h for which the trivalent 
shell has lower energy than both the divalent and tetrava- 
lent shell. Experimentally, however, we do observe diva¬ 
lent, trivalent and tetravalent shells. Furthermore, they 
appear with similar likelyhood with a tendency towards 
formation of the trivalent shell, as shown in the histogram 
of Fig. [^for shells with h/R 0.1 made using microflu¬ 
idics. The presence of all shell types likely result from 



FIG. 6 : Elastic energy as a function of shell thickness 
for divalent (red), trivalent (green) and tetravalent 
(blue) defect configuration for Rja = 10^ and Ec = 0. 
Either the divalent or tetravalent configuration, but not 
the trivalent configuration, is lowest in energy. 
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EIG. 7: Goexistence of divalent, trivalent and 
tetravalent configurations in shells with hjR ^.1. The 

histogram shows the different populations in a sample 
of 60 shells. 


their similar energies; this makes all three of them acces¬ 
sible while they are being generated. The lack of shell 
type transitions and hence the observed lifelong stability 
of a certain shell type irrespective of whether it corre¬ 
sponds or not to the ground state is likely due to en¬ 
ergy barrier between shells with different defect number. 
To ascertain this, we calculate the corresponding energy 
barriers. Eor going from the trivalent to the tetravalent 
shell, we simply have to undo the escape in the third di¬ 
mension. The associated barrier thus corresponds to the 
difference between ^^=3 and Ez= 3 - 


AE3^4 = irKh 




-4.19 


Ec 

TvKh 


(16) 
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The energy barrier for going from the divalent to the 
tetravalent shell lies in overcoming the repulsion between 
the two + 1/2 defects and can thus be estimated by the 
difference between E '^^2 ^z= 3 ’ 


AJ53..=,A-;.(ll„g(5)-0.I5-(17) 

Since K « for a thin shell with h = 1/im, 

these barriers are four orders of magnitude larger than 
the thermal energy scale (at room temperature), ksT^ 
where T denotes temperature and is Boltzmann’s con¬ 
stant, thus providing stability of any of the shell types 
after made, explaining why valence transitions are not 
observed experimentally. 

Finally, to understand the tendency towards the triva- 
lent shell, we recall that while the shells are generated, 
the inner droplet is usually never at the center of the 
outer liquid crystal droplet. This implies that the shells 
are heterogeneous in thickness, with a thinner and a 
thicker part. Since defect nucleation happens at random 
locations inside the shell and for the thin part h < ^ 23 , 
where E 4 < Es < E 2 , we might be providing conditions 
for formation of two + 1/2 lines in this region. Simi¬ 
larly, since in the thick part of the shell, h > ^ 34 , where 
< ^3 < ^ 4 , we might be providing conditions for for¬ 
mation of a pair of +1 boojums in this part of the shell. 
The tendency observed experimentally towards the triva- 
lent shell could then relie on the thickness heterogeneity 
of the shells as they are made in our microfluidic devices. 
Once the shells are made, their thickness and thickness 
inhomogeneity reach their equilibrium configuration and 
the defects locate at the top, where the elastic free energy 
of the nematic liquid crystal is minimum. 


to rewrite E in eq. © as an expansion to quadratic order 
in q\ 


F = F {el 4>1 + \Y^ Mijmj + O (q^) . (21) 

ij 


The 6 X 6 matrix M can thus be found by 


Mij — 


d‘^E 

dqidqj 


q-=qj=0 


( 22 ) 


This calculation is performed without loss of generality 
upon choosing the ground state defect locations to be on 
the equator, i.e. 0^ = 7 r/ 2 . We diagonalise this matrix: 


M = UDU^. (23) 


The matrix D has the following eigenvalues on the diag¬ 
onal: 


{AJ 


nk 


{0,0,0,15,17,18} 


(24) 


The columns of the matrix U are the corresponding or¬ 
thonormal eigenvectors, {u^}, and is the transpose of 
U. The eigenvectors belonging to the three zero eigen¬ 
values represent rigid body rotations. The other eigen¬ 
vectors are 


U 4 


0 

0 

_]_ 

+2 

V 72 / 


U5 



3 

Vm 

3 

VM 

0 

0 


V 0 / 


U6 


/ 0 \ 


0 

0 



(25) 


IV. BOND FIDELITY 

In this section, we will consider the fidelity of the three 
bonds by considering its robustness against thermal fluc¬ 
tuations. We will expand the energy around the equilib¬ 
rium values for the zenith and azimuthal angles, { 6 >/, 0 ^}. 
We parametrise the departures from the equilibrium an¬ 
gles with a 2 Z-component vector q, whose first three 
components are the deviations along the lines of longi¬ 
tude of the sphere and whose final three components are 
the deviations along the lines of latitude of the sphere. 
We thus have 

qi = S 0 u (18) 

qs+i = S(l)ism0i. (19) 

Again we employ the law of cosine on the sphere 


The fourth and sixth eigenvalues also correspond to de¬ 
formations that keep the defects located at a great cir¬ 
cle. The fourth one corresponds to a displacement of the 
charge one-half defects such that their distance to the 
charge one defect grows or shrinks in equal manner and 
hence preserves the isosceles shape of the triangle (Fig. 
§ 1 ). The sixth eigenvalue corresponds to a mode defor- 



FIG. 8 : Schematics of the three non-trivial eigenmodes 
corresponding to (a) U4, (b) U5 and (c) uq. The defects 
(represented by dots) continue to lie on a great circle in 
(a) and (c), but not in (b). The defects continue to 
form an isosceles triangle in (a) and (b), but not in (c). 


cosPij = cos ( 6 >/ + qi) cos [Oj + qj) + sin ( 6 >/ + qi) sin ( 6 >° + 


qs+i 
sin 6 >? 


g3+j 

sinO^ 


( 20 ) 


X cos 


mation that does not posses this property, thus breaking 
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the symmetry of reflection of the bisector of the distinct 
angle (Fig. [^). The mode of deformation correspond¬ 
ing to the fifth eigenvalue, however, retains the isosceles 
shape of the triangle, but shrinks the size of the trian¬ 
gle as the defects do not lie on a great circle anymore 
(Fig. [^). We change the basis from qi to Wi^ which is 
the departure from the trivalent ground state in the Tth 
eigendirection: 

Qi = UijWj. (26) 

This basis transformation yields to quadratic order in wf. 

F = F (^0^, -X 4 w‘lX qWq. (27) 

By equipartition, each term contributes ^ksT. The 
eigenvalues corresponding to the trivalent modes of defor¬ 
mation are equal or larger than the the tetravalent ones 
(which are and ^nk [2 [5]): the trivalent ground 
state is thus somewhat better protected against thermal 
fluctuations. 


boojurns is observed experimentally besides the bipolar 
and regular tetrahedral configuration. The repulsive in¬ 
terdefect interaction pushes the defects to lie on a great 
circle. The strength of the interaction depends on the 
charges of the defects. Consequently, the defects are lo¬ 
cated at the vertices of a isosceles triangle rather than a 
equilateral triangle, in contrast to the tetravalent ground 
state in which the defects are equidistant. In the energet¬ 
ically most favourable trivalent configuration, we obtain 
for the central angles ^12 = ^13 = 0.737r, (323 = 0.547r and 
for the angles in the (flat) isosceles triangle ai = 48° and 
0:2 = 0:3 = 66°. These values are in agreement with ex¬ 
perimental values. Estimations of the elastic energy show 
that there is no shell thickness for which the trivalent 
ground state is lower than both the tetravalent and diva¬ 
lent ground state. However, there are energy barriers to 
provide stability for the trivalent state once it is created. 
We note that our calculations do not include thickness 
heterogeneity. Including this effect is the next step to¬ 
wards our continued understanding of trivalent nematic 
shells. 


V. CONCLUSION 

In a spherical nematic shell of finite thickness a stable 
defect structure with two s = 1/2 lines and one pair of 
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